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ABSTRACT 

Early-type dwarfs (dEs) are by far the most abundant galaxy population in nearby 
clusters. Whether these objects are primordial, or the recent end-products of the differ- 
ent physical mechanisms that can transform galaxies once they enter these high-density 
environments, is still a matter of debate. Here we present a novel approach to test 
these scenarios by comparing the properties of the globular cluster systems (GCSs) 
of Virgo dEs and their potential progenitors with simple predictions from gravita- 
tional and hydrodynamical interaction models. We show that low-mass (M* < 2 x 10 s 
Mq) dEs have GCSs consistent with being the descendants of gas-stripped late- type 
dwarfs. On the other hand, higher mass dEs have properties -including the high mass 
specific frequencies of their GCSs and their concentrated spatial distribution within 
Virgo- incompatible with a recent, environmentally-driven evolution. They mostly 
comprise nucleated systems, but also dEs with recent star formation and/or disc fea- 
tures. Bright, nucleated dEs appear to be a population that has long resided within the 
cluster potential well, but have surprisingly managed to retain very rich and spatially 
extended GCSs - possibly an indication of high total masses. Our analysis does not 
favour violent evolutionary mechanisms that result in significant stellar mass losses, 
but more gentle processes involving gas removal by a combination of internal and ex- 
ternal factors, and highlights the relevant role of initial conditions. Additionally, we 
briefly comment on the origin of luminous cluster SOs. 

Key words: galaxies: clusters: individual: Virgo - galaxies: formation - galaxies: 
evolution - galaxies: dwarf - galaxies: interactions - galaxies: star clusters: general 



1 INTRODUCTION 

Early-type cluster dwarfs^ constitute a longstanding puzzle. 
Being the mos t abundant galaxy p opulation in high-density 
environments (|Binggeli et alj |1985). it is not yet fully under- 
stood whether they are old systems whose current properties 
were determined at the early stages of cluster assembly, or 
if they have a late origin as the end-products of the physical 
mechanisms that operate in virialised clusters. The scenario 
is indeed complex, as it is well known that dEs are not sim- 
ple objects with uniform properties, but rather encompass 
a variety of systems with different stellar structure and age 
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by this definition we refer to 



any low-mass, quiescent cluster galaxy, thus including dwarf el- 
lipticals, spheroidals and lcnticulars. Throughou t this paper they 
will all be termed dEs for simplicity - but sec Kormcnd v et all 
for a thorough discussion of the similarities and differences 
between these specimens. 



(e.g., iFerrarese et~afl l200rj : iBoselli et all 120081 ; iLiskerl 120091 

and references therein). It is therefore most likely that sev- 
eral formation mechanisms need to be invoked to explain 
the origin of these objects. 

These mechanisms can be broadly divided into inter- 
nal and external ones. Among the first class, strong feed- 
back effects due to supernovae and intense star for mation 
is arguably the best candidate ( Dek el fc Silkl Il986j ) , as it 
can lead to the removal of low angular momentum material 
-be it dark matter (DM) or the interstellar medium (ISM)- 
, potentially resulting in a thick system that is not DM- 
dominated in the central regions b ut can still be rotationally 
supported (|Governato et al1l2010h . The distinct location of 
dEs on the Fun damental Plane app ears to be consistent with 
this scenario l|de Riicke et al.ll2005l) - although environmen- 
tal eff ects can not be ruled out ( Aguerri fc Gonzalez-Garcial 
l2009h . 

The second class includes a plethora of physical 
mechanisms that can shape galaxies in clusters (see 
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Bos elli fc Gavazzil l2006 for an excellent review). These ex- 
ternal processes can be again subdivided in two groups: 
hydrodynamical interactions between the hot intracluster 
medium and the gala xy's ISM -e.g., r am-pressure strippi ng 
l|Gunn fc Gottlll972T ) and starvation l|Larson et alJll980h - 
and gravitational interactions with the cluster potential and 
between galaxies th emselves -tidal str ipping (|Merritdll984l ) 
and tidal shocking |Moore et al.|[l9 96). In reality, all these 
mechanisms contribute to some extent to the evolution of 
cluster dwar fs, in a sort o f tidal stirring process like the one 
proposed bv lMaver et al.l (|200lh . 

Environmental effects are supported by a significant 
number of observational evidences, including the fact that a 
non-neg ligible fraction of dEs appear to be rotationally sup- 
ported JPedraz et al.ll2Q02l ; iGeha et al.l 12003 ; Ivan Zee et al.l 
12004 iToloba et al.ll201ll ); the presence of disc features (spi- 
ral arms and/or bars) in a subpopulation of cluster dEs 
(Bar azza et al. 20021; iGraham fc Guzma n 2003; Liskcr ~t al.l 



120061 ; iMendez-Abreu et al.ll2010l ; Ijanz et alj|2012h : the ex- 
istence of both a bulge and disc com ponents in th e light 
profiles of some early-ty pe dwarfs (|Aguerri et alj 120051 ; 
iKormendv fc Benderll2012i); or the strong transit ion towards 
older ag es ilHaines et al.ll2006l;ISmith et alj|2009h and redder 
colours (jSanchez- Janssen et alj 20081 ) of dwarfs within vit- 
alised cluster regions. Furthermore, the signatures of a late 
infall from cluster dw arfs dynamics (|Drinkwater et aT1l200ll ; 
IConselice et al.ll200ll ) and the realisation that the low-mass 
end of the red sequence pro gressively builds up at lower red- 
shifts (|De Lucia et al.|[2007l ) strongly suggest that this trans- 
formation is a relatively recent phenomenon that took place 
over the last ~ 6 Gyr. Unfortunately, all these mechanisms 
can be efficient in multiple, sometimes coincidental cluster 
regions, and their timescales are also different. To observa- 
tionally discriminate the predominant process(es) becomes 
a difficult task, requiring the use of many observables that 
can help break these degeneracies. 

Globular clusters (GCs) are ancient systems found in al- 
most all galaxy types along the Hubble s equence, irrespec- 
tive of their luminosity or gas content (jBrodie fc Straderl 
120061 . BS06 hereafter). Thus, the ubiquity of these stel- 
lar clusters have made them powerful tools for the study 
of star formation and mass assembly processes in nearby 
galaxies, as they hold a fossil record of their host galaxy 
evolutionary history. This fact has prompted the use 
of globular cluster systems (GCSs) as sensitive tracers 
of formation mechanisms and, to a lesser extent, envi- 
ronm ental effects on different cluste r galaxy popu l ations 
(e.g.. lAragon-Salamanca et al.1 |2006l: IMiller fc Lotj 120071 ; 



ICoenda et all l2009h . iPeng et all (|2008l ) provide the most 
comprehensive analysis to date of the GCSs in Virgo dE 
galaxies, including an investigation of the roles of stellar 
mass and environment. They show that dEs in the inner 
cluster regions exhibit higher GC specific frequencies than 
those in the outskirts, but are indistinctive in terms of their 
structural properties or their stellar populations colours - 
a behaviour they explain within a framework of an early 
origin for dEs and a biased GC formation towards dense 
environments. 

In this contribution we extend their approach by quanti- 
tatively comparing the properties of the globular cluster sys- 
tems of Virgo dEs and their potential progenitors with sim- 
ple predictions from the gravitational and hydrodynamical 



interaction models that are traditionally thought to shape 
early-type cluster dwarfs. The Paper is organised as follows. 
In Section 2 we present a brief overview of the data samples 
compiled from the literature that are used during the sub- 
sequent analysis. Section 3 presents predictions of the joint 
evolution of galaxies and their GCSs from the most relevant 
interaction models in clusters. Sections 4 is devoted to the 
constraints derived from observations, which are further dis- 
cussed in Section 5, together with future prospects that we 
consider can significantly improve our knowledge of dwarf 
galaxy formation and evolution. 



2 DATA COMPILATION 

In order to carry out this study, we have compiled data for 
the GCSs of cluster and field galaxies from several sources 
in the lite rature. Most ob servational data come from the 
ACSVCS (|C6te et alj|2004 ). an unprecedented HST survey 
of 100 early-type galaxies in the Virgo cluster spanning a 
wide magnitude range ( — 23 < Mv < —15.5). Taking ad- 
vantage of the HST's unique resolution, they have carried 
out a thorough analysis of the GCSs in their sample, in- 
cluding st udies of GC colo u r dist ributions and s pecific fre- 
quenc ies (|Peng et alj|2006t I2008T ). Furthermore. iMei et ail 
l|2007l ) derived surface brightness fluctuations distances for 
the majority of the sample, allowing for a projection- free 
investigation of environmental effects as a function of posi- 
tion within Virgo. Stellar masses, derived from colour fits to 
m odel single stellar populations, are also used as presented 
in lPeng et all {2008, P08 hereafter), and galaxies having M* 
< 5 x 10 9 M are considered to be dwarfs. 

In order to extend the analysis further into the low-mass 
regime, we have complemented the previous d ataset with the 
faint (-18 < M v < -12.5) Virgo dEs from IMiller fc Lotzl 
i|2007l . ML07 hereafter). A mass-to-light ratio in the U-band 
of M*/W = 2 is used when computing their stellar masses, 
consistent with those of the ACSVCS. 

The properties of the GCSs of potential progenitors 
come from several so urces. At the bright end we use the 
recent compilation of ISpitler et all l|2008l , S08 hereafter), 
that includes data for several M* > 4 x 10 10 M Q ellipti- 
cal and disc galaxies. This dataset is particularly well suited 
for our analysis, as it only contains galaxies having total 
GC number uncertainties lower than 40%, a large spatial 
coverage (more than 50% of the GCS radial extent), avail- 
able information on metallicity subpopulations proportions, 
and self-consistently derived stellar masses. We further com- 
pleme nt this dataset with the late-type galaxy compila- 
tion of lChandar et all l|2004h . using the morphological type- 
dependent M t /Lv ratios derived by S08. Given that the 
majority of galaxies in the sample lack information regard- 
ing subpopulation fractions, we assume they all consist of 
metal-poor GCs. These values are therefore much more un- 
certain than the S08 ones (and, probably, upper limits), but 
we nonetheless include them in our analysis for comparison 
purposes. 

At the faint end we use data from the recent work of 
iGeorgiev et all (|2010l . G10 hereafter), consisting of low-mass 
(M* < 3 x 10 9 M Q ) galaxies in the Local Universe (D < 
10 Mpc). From this work we only select late-type (T 5) 
dwarfs, and use stellar masses as tabulated by the authors. 
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A dditionally, we als o include the seven Virgo dlrrs studied 
bv lSeth et all |2004 S04), for which IVL/W = 1 is assumed 
- consistent with both the G10 and P08 values. 

We note that these samples comprise the most homoge- 
neous and robust compilations of GCSs for their respective 
galaxy types. This is especially true in the case of dwarf 
galaxies, where all four studies benefit from the unrivaled 
resolution and sensitivity of the HST. 



3 ENVIRONMENTAL EFFECTS ON 
GALAXIES AND THEIR GLOBULAR 
CLUSTER SYSTEMS: MODEL 
PREDICTIONS 

Globular clusters, with their predominantly old agesQ, have 
been privileged witnesses of a significant fraction of their 
host galaxy evolutionary history. This, together with their 
intrinsic compact sizes and extended galactocentric distri- 
butions, makes them ideal probes of interaction processes in 
galaxy clusters. For instance, they are in principle insensi- 
tive to all the hydrodynamical mechanisms that can alter 
the properties of a galaxy's ISM once it enters the cluster 
environment. On the other hand, gravitational interactions 
can modify the stellar content of galaxies, and this should 
be reflected on their GCSs as well. Here we make predic- 
tions of the joint evolution of galaxies and their GCSs using 
basic theoretical considerations of the most relevant inter- 
action models in clusters, i.e., gravitational harassment and 
gas stripping. We focus on relatively recent interactions (< 6 
Gyr), and implicitly assume that none of these processes lead 
to GC formation. 



3.1 Gravitational Harassment 

One of the most popular evolutionary paths for clus- 
ter dEs is their transformation from late-type galax- 
ies through multiple, fast gravitationa l encounters withi n 
the cluster potential, i.e., harassment (|Moore et al.l fl99"rjh . 
This tidally-induced mechanism is supposed to be rather 
efficient at stripping gas and stars from both rela- 
tively massive, lo w surface brightness galaxies (LSBs; 
M oorc et all Il999l ) a nd lower-mass dwarf irre gular (dirr) 
and Magellanic- types (|Mastropietro et al.ll2005l ) . The result- 
ing remnants closely resemble o bserved cluster dEs (e.g., 
lAguerri fc Gonzalez-Garci'al 120091 ). 

Gravitational interactions in clusters like Virgo are suf- 
ficiently rapid t hat they can be treated under the impulse 
approximation (Spitzcr 1958). These fast encounters result 
in a net increase of the kinetic ene rgy per unit mass of th e 
target's particles in the form (e.g.. lAguilar fc Whitel ll985): 



(AE/m) = G 2 M% v~ 2 b~ 4 r 2 f(P, A), 



(1) 



2 The simulations of IBeaslev et all ll2002t ) predict mean ages of 
12 and 9 Gyr for the metal-poor and metal-rich subpopulations, 
respectively - though the latter are characterised by a much 
broader age distribution. This is consistent with observations of 
both Galactic and extragalactic GCs (see BS06 and references 
therein). In the particular case of the Virgo cluster, old GC ages 
(>10 Gyr) have been inferred from spectroscopic observations for 
galaxies spanni ng two orders of m agnitude in stellar mass, from 
the giants M87 dCohen et al.lll998l) and M 49 llPuzia et al.lll999l) . 
to the dE VCC 1087 [IBeaslev et alJlioOrJ) . 



where G is the gravitational constant, Mp is the perturber's 
mass, v is the encounter velocity, b is the impact param- 
eter, and r is the distance of the affected particle to the 
target's centre. This energy change is further modulated by 
the f(P, A) term, that includes a correction for both the 
extended nature of the perturber, f{P), and for the move- 
ment of the stars within the target - the adiabatic correc- 
tion. The latter can take the form f(A) — (1 + x 2 )~"' , where 
x — lot is the product of the orbital frequency of parti- 
cles in the target (oj) and the effective duration of the en- 
counter (r), and 1.5 < 7 < 3 depending on the duration 
of the sho ck with respect to t he target's half-mass dynam- 
ical time (|Gnedin et al.| [l999). Equation[T] implies that the 
energy gain is maximal for particles at large galactocentric 
distances with small orbital frequencies, while it is strongly 
diminished if they are concentrated or if tidal forces act on a 
time scale much larger than their orbital period. The larger 
the kinetic energy acquired, the less bound particles are. As 
a result, the strong dependence on galactocentric distance 
(ocr 2 ) implies that in the impulsive regime mass loss pref- 
erentially occurs from the outside-in. 

The efficiency of mass stripping due to harassment de- 
pends on the collisional frequency, and on the s trength of in- 
dividu al collisions and of the cluster tidal field. iMoore et al.l 
( 1998) show that the competing effects of these terms re- 
sult in harassment being a rather effecti ve mechanism in a 
wide range of clustercentric radii (see also lGnedinll2003h ■ but 
more so closer to the cluster centre, where both the cluster 
tidal field and the collision rate are expected to be maxi- 
mal. Taking all this into account, we can make the following 
predictions for the evolution of GCs and the stellar body in 
harassed galaxies: 

i) The detailed stripping history of both components will 
be determined by their phase-space (spatial and orbital) 
distribution. If GCs have a more extended spatial distri- 
bution and lower rotational suppo rt than stars, they should 
be stripped more efficiently (e.g., Be kki et a L 20031 ) - and 
viceversa. 

ii) The effects of mass stripping should present a cluster- 
centric dependence, in the sense that objects that spend 
longer times close to the centre are affected the most (e.g., 
iMastropietro et alJliooBI : ISmith et al]|2010l ). 



3.2 Gas Stripping 

The alternative and perhaps most obvious scenario for 
the origin of dEs is gas removal from a population of in- 
falling, field dlrrs. Gas sweeping can occur through several 
processes that act on differen t timescales, i ncluding ram - 
pressure dGunn fc Gottl Il972l) and visc ous i|Nulsenl Il98^ ) 



strippin g, starvation jLarson et al.lll980h and thermal evap- 
oration (|Cowie fc Songailal 1977 ). They all can halt star for- 
mation and -especially ram-pressure- create low-mass, pas- 
sive galaxies structurall y and chemically c ompatible with 
those observed in Virgo l|Boselli et al.l [20081 ) . 

Whatever the specific mechanism, it is important to 
note that their effect o n the already-exist ing stellar body 
is essentially negligible (iQu ilis et alj 12001: i , and the same 
is true for the GCS. Therefore, if gas-stripping mechanisms 
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play a predominant role we expect to see no difference what- 
soever between the GCSs of potential progenitors and rem- 
nants. One important caveat is that this is not exactly true 
for their stellar mass content. Indeed, while star formation 
in field dlrrs has continuously increased M*, in the gas- 
stripping plus star formation shutoff scenario the current 
masses of dEs should essentially reflect those of their pro- 
genitors upon cluster arrival. This stellar mass difference will 
generally depend on the detailed star formation history of 
the progenitor, and on the efficiency and timescale of the 
stripping process - resulting in a 'progenitor bias' that has 
to be taken into account. 



non-monotonic relation between TN,blue and galaxy stellar 
mass is clearly depicted in this plot: intermediate-luminosity 
galaxies exhibit rather uniform, low values (note the loga- 
rithmic scale) , while there is a clear trend for both low-mass 
and high-mass galaxies to have increasing specific frequen- 
cies. This already hints at the possibility that a transition 
from the intermediate-mass to the low-mass regime most 
likely requires an increase in Tn, Hue- We develop on this 
scenario by providing a more rigorous comparison between 
the GCSs of cluster dEs and their potential progenitors. It 
is important to note that evolutionary paths that decrease 
(increase) the number of GCs of a galaxy necessarily have 
shallower (steeper) slopes than the constant N gc diagonal 
lines in Fig.[T] 



4 CONSTRAINTS FROM OBSERVATIONS 

4.1 GC Mass Specific Frequencies 

One of the most constraining observables provided by GCSs 
is their specific frequency, i.e., the number of GCs nor- 
malised by galaxy luminosity or s tellar mass. If GCs indeed 
form with a universal efficiency (|McLaughlinl 1 19991 ) . varia- 
tions of specific frequency among galaxy types might reflect 
a different evolutionary history. For example, if dEs are the 
descendants of transformed galaxies, the properties of both 
GCSs must be somehow related - and if several evolutionary 
scenarios lead to the formation of these objects we expect 
their distinct imprint to be recognisable on their GCSs. 

Here we explore the connection between dEs 
and their potential progenitors by means of the GC 
mass specific f requen cy, Tn = 7V gc /(M*/10 9 M@). 
IZepf fe Ashman! J 19931) modified the c lassic Sn defini- 
tion of lHarris fc van den Berghl (|198lh -replacing the 
galaxy luminosity by its stellar mass-, thus allowing for 
a more straightforward comparison between galaxies with 
different M*/L ratios. Moreover, this variation turns out to 
be a better discriminator between interaction types than 
Sn- From the parameters involved in these definitions, it is 
clear that both Tn and Sn are expected to vary during a 
gravitational encounter involving any kind of stellar mass 
loss (stars or GCs) . Neither the stellar mass of a galaxy nor 
its number of GCs are modified during a gas stripping event 
- thus leaving Tn unchanged. But the lack of gas reservoirs 
for future star formation unavoidable leads to luminosity 
fading -and Sn varying- due to stellar populations aging. 
In other words, Sn variations are expected due to both 
gravitational and hydrodynamical processes, while Tn is 
only sensitive to the former. 

During this discussion we further focus on the mass spe- 
cific frequency of metal-poor (blue) GCs. This subpopula- 
tion, unlike the metal-rich (red) one, is unambiguously old 
and appears to be ubiquitous across galaxies of all masses 
and morphological types - therefore allowing for a direct 
comparison between dEs and their potential progenitors. 
Furthermore, they are traditionally associated to the halo 
of L > L» galaxies due to their extended spatial distribu- 
tion and kinematics (BS06), making them excellent tracers 
of gravitational interactions (cf. Sect. l3.l| l. 

In Fig.[T]we present the distribution of Tjv,w«e for metal- 
poor clusters as a function of M* for Es, SOs, Sps and dEs in 
our compilation. Diagonal lines correspond to the location 
of objects having from one to 10 4 blue GCs. The well-known 



4.1.1 LSBs vs dEs 

We start by comparing the GC mass specific frequency of 
LSBs and dEs with a sim ple scaling of the harassment mod- 
els of lMoore et~aT1 1| 19991 . M99 hereafter). These authors fol- 
low the dynamical evolution of a M* ~ 10 10 Mq LSB within 
a massive galaxy cluster. Given that GCSs were not included 
in the simulations of M99, we estimate their evolution as fol- 
lows. 

The GCSs of luminous disc galaxies are characterised 
by having rather homogeneous properties. Their specific fre- 
quen c y appears to be alm ost universal (|Goudfrooii et all 
120031 ; Ichandar et al.l |2004| ; see also Fig.[l]), and they have 
a more extended spatial distribution than disc stars - 
especially the metal-poor population (see BS06). Interest- 
ingly, these features appear to be independent of galaxy sur- 
face brightness, as indicated by the an alysis of GCSs in five 
galaxies spanning a wide range of fio (|Villegas et al.l [20081 . 
V08 hereafter). We stick to these results and assume that 
the properties of the GCSs of intermediate-luminosity LSBs 
are similar to those of high surface brightness (HSBs) disc 
galaxies like M31 or the Milky Way.0 Fig.[T] illustrates the 
narrow range of frequencies cha racteristic of disc ga laxies, all 
having 1 < Tn^Iub < 4 (see also lRhode et al.ll2007r ). In order 
to match the general trend in Fig. [T] we assign our 10 10 M Q 
model LSB a fiducial mass specific frequency TN,biuc ~ 5 
(large solid star) - or, equivalently, an spatially extended 
population of N gc ,blue = 50. This is a generous upper limit 
to the metal-poor GC specific frequencies observed in disc 
galaxies. M31 and M87, with TNfilue ~ 2 and 35, respec- 
tively, are labeled for reference. 

We noted in Section l3~T1 that the detailed stripping his- 
tories of stars and GCs will critically depend on their spatial 
and kinematical properties, which are poorly constrained 
for disc galaxies outside of the Local Group. Fortunately, 
we can combine the expectations from the impulse approx- 
imation and results from numerical simulations to provide 
an order-of-magnitude estimate of their evolution. Observa- 
tions indicate that the metal-poor GC population in L* disc 
galaxies has a more extended spatial distribution than disc 
stars. Additionally, they most likely are less extended than 



3 We stress that our assumption only concerns the similarity of 
the GCS and not of the galaxy as a whole. M99 show that HSBs 
barely lose stellar mass through harassment and therefore can not 
evolve into dEs - but they might do into SOs (see AppendixlATl . 
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Figure 1. Mass specific frequency of blue (metal-poor) GCs as a function of stellar mass for E, SO, Sp and dE galaxies from the literature. 
Galaxies having no GCs have been assigned a Tn blue = 0.12, and diagonal lines indicate the loci of galaxies having from one to 10 4 blue 
GCs. The solid star shows the initial Tff j,i ue of our fidu c ial LS B model, while the shaded region indicates its potential location according 
to our scaling of the harassment model of iMoore et al.l lll999l) . The blue curve indicates its expected evolution if mass stripping occurs 
from the outside-in (see text for details). The only early- type dwarfs compatible with being harassment remnants are those showing low 
Tjv blue values. LSB galaxies would need to lose more disc stars than GCs in order to match the high specific frequencies observed in 
most dEs. 



particles of the DM halo. These two constraints to the metal- 
poor GCS spatial distribution translate into the expectation 
that GCs should be more easily stripped than stars, but less 
than DM. Moreover, this probably holds even at a fixed 
radius, as kinematical studies of the best studied disc galax- 
ies -M31 and the MW, see BS06- indicate that GCs move 
on more (less) eccentric orbits with less (more) rotational 
support than disc stars (DM). Thus, we generically predict 
that the fraction of stripped GCs will always be intermedi- 
ate to that of the DM halo and the bulk of the stellar disc 
(fdm ^ fgc ^ /disc)- In order to provide a more quantitative 
comparison with observations, we use the results from the 
harassment simulations of M99 (their table 1), where it is 
shown that up to 90% of disc stars and 95% of DM mass 
can be stripped - i.e., < fdisc < 0.9 and < fdm < 0.95. 

The solid star in Fig.[T] shows the initial specific fre- 
quency of the model LSB, while the shaded region indi- 
cates its potential location according to our evolutionary 
predictions. The horizontal upper limit corresponds to the 
case where the LSB loses the same fractions of GCs and 
disc mass, hence T^ t uue remains constant. The rightmost 
and leftmost vertical boundaries indicate the extreme cases 
where the stellar mass losses are 0% and 90% of the origi- 



nal value, respectively - but the fraction of stripped GCs is 
always larger and therefore Tisr,biue decreases. The lower di- 
agonal boundary is the limiting case where the galaxy is left 
with one GC only, while the stellar mass loss ranges between 
the two chosen extremes. Interestingly, the low-mass, high- 
TN,uue region of the diagram is not accessible by this type 
of interaction, implying that such dEs are difficult to under- 
stand within the harassment scenario. The only way an LSB 
remnant can increase its TM,Mue is by preferentially losing 
more disc stars than GCs, an unlikely evolution according to 
the impulse approximation. Alternatively, LSBs must have 
very different GCS properties (abnormally high frequencies 
or a very concentrated spatial distribution) than the bulk 
of the disc galaxy population - but this is not supported by 
observations. 

The blue curve in Fig.[T] provides a more specific pre- 
diction of the evolution of harassed LSB galaxies. We ex- 
amine the limiting case where stripping occurs strictly from 
the outside-in, i.e., without taking the orbital properties of 
stars and GCs into account. The model consists of a scaled 
version of the Milky Way. We assume that the metal-poor 
clusters have the same spatial distribution as in the MW - 
a standard practice in studies of extragalactic GCSs in disc 
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galaxies (e.g., iGoudfrooii et all 120031 ; IChandar et "all 12004 
V08). We use the Galactic GCS data compiled by Bica et all 
(2006), and select metal-poor clusters as those with [Fe/H] 
s£ -1. The LSB model has a M* = 10 10 M Q disc with the 
same exponential scale length as the MW (h = 2.5 kpc) - 
and therefore, a factor 10 lower central surface mass den- 
sity. At any given radius we estimate the fractions of en- 
closed GCs and disc mass, and the (properly normalised) 
ratio provides the evolution of the mass specific frequency. 

As expected, Tn,blue rapidly decreases at first due to 
the outermost GCs being stripped without significant stel- 
lar mass loss taking place. Stripping then continues with 
more or less equal efficiency for both components, result- 
ing in a low-mass, low-Tjv,w U e remnant. This is admittedly 
a simplistic approximation, as we assume that the spatial 
distribution of stars and GCs are not altered during the evo- 
lution, but nonetheless shows that it is extremely difficult to 
increase Tmmub if stripping occurs from the outside-in. It is 
of course possible that, as mass loss proceeds, GCs in the in- 
ner regions start to become more resilient to stripping than 
stars. What Fig.[T] shows is that at this point the number of 
surviving GCs is already too low (N gc ^iue J; 10) to match 
the rich GCSs observed in most massive dEs. And, given 
that metal-poor GCSs tend to be less rotationally supported 
than disc stars (BS06), we can further expect this Tn,blue 
evolution to be an upper limit to the actual process. 

We finally note that, most likely, the GCSs of Virgo 
dEs have also been affected by the cluster tidal field dur- 
ing their recent orbital evolution, resultin g in some level o f 
stripping and/or enhanced orbital decay l|Oh fc Li"nl [20001 
If their GCSs were richer in the past, the mass specific fre- 
quencies currently observed can be thought of as lower lim- 
its - thus making the disagreement with disc galaxies even 
stronger. 



4- 1.2 dE subpopulations , dlrrs, and the importance of 
environment 

Recent observational work indicate that, on average, early- 
type dwarfs have much richer GCSs than late- types (e.g., 
ML07, G10). This is illustrated in Fig.[2l a zoomed-in version 
of Fig.[T]where we show the mass specific frequency of metal- 
poor GCs for Virgo dEs (ML07, P08) and dlrrs (S04), and 
field dlrrs from G10. There are several features in this Figure 
worth discussing. 

First, the richest GCSs among dEs are always found at 
the high-mass end: essentially all galaxies having more than 
ten blue GCs are more massive than M* ~ 2 x 10 8 . This re- 
sults in a significant GC abundance difference with respect 
to field dlrrs, where only 2/33 reach that number. Remark- 
ably, dlrrs in the Virgo cluster appear to be richer than their 
counterparts in low-density environments, and they nicely 
match the high specific frequencies of dEs. In the narrow 
5 x 10 7 < M^/Mq < 5 x 10 s mass range, field dlrrs have a 
median number of N gCt ti U e = 3 metal-poor GCs, compared 
to Ngc,blue — 13 for Virgo dlrrs (see Fig.|2J. This difference 
was first uncovered by S04, and argues in favour of a funda- 
mental difference between field and cluster late-type dwarfs 
in general. For instance, NGC 1427A -the dlrr galaxy hav- 
ing the richest GCS in Fig. [5] from the G10 sample- actually 
is a member of the Fornax cluster. 

We recall that a direct comparison between early- and 



Table 1. Statistical properties of the metal-poor GCSs in the dif- 
ferent M* > 2 X 10 s Mq dwarf types in Fig.[2] Each non-diagonal 
element in the upper table indicates the confidence level, accord- 
ing to a Kolmogorov-Smirnov test, with which we can rule out 
the null hypothesis that the two Tjv blue distributions are drawn 
from the same parent population. The diagonal elements show 
the fraction of galaxies of each type that have ten or more blue 
GCs. In the lower table we show the corresponding lower quar- 
tiles, medians and upper quartiles of the Tjy blue distributions. 



/(> 10 GCs) 


dlrr 


dE(di,bc) 


dE(nN) 


dE(N) 


dlrr 


12 


78.8 


97.5 


99.8 


dE(di,bc) 


78.8 


46 


63.4 


92.3 


dE(nN) 


97.5 


63.4 


50 


80.3 


dE(N) 


99.8 


92.3 


80.3 


78 


Ql(TN,blue) 


3 


5 


7 


11 


Qz(Tn ,blue) 


5 


8 


12 


19 


Q3(TN,blue) 


10 


18 


18 


33 



late-type dwarfs is hampered by the different star forma- 
tion histories (i.e., stellar mass growth) of these populations 
(cf. Sect. l3.2| l. To take this potential 'progenitor bias' into 
account, the arrow in Fig. [2] indicates the effect of M* dou- 
bling had the dE progenitors not had their star formation 
quenched. In this case galaxies simply would move along 
constant N gc paths, and this does not change the fact that 
most dEs contain way too many GCs compared to field dlrrs. 

We have so far not distinguished between the differ- 
ent subpopu lations in which V irgo dEs can be divided. As 
discussed bv lLisker et all l|2007l ). there are at least four sub- 
classes of dEs with distinct morphological and clustering 
properties: dE(bc), with central star formation; dE(di) ex- 
hibit disc features; dE(nN) have weak or no nuclei; and 
the classic nucleated dE(N). Put in this order, they can be 
naively thought of as forming a temporal sequence of cluster 
membership, from the most recent newcomers to the oldest 
cluster inhabitants. It is therefore pertinent to investigate 
whether these different properties also extend to their GCSs, 

as testament to their evolutionary histories. 

In Fig. [2] we adopt the classification of iLisker et all 

|2007h for each dE, but group dE(bc) and dE(di) together 
for the sake of statistics. Even though the ACSVCS imag- 
ing has demonstrated the presence of nuclei in galaxies that 
were thought no t to have one from ground-based studies 
(|C6te et al.ll2006l ). here we maintain the VCC nuclear clas- 
sification for homogeneity reasons - ML07 did no t re-cl assify 
their galaxies. We therefore stick to lLisker et all l|2007h and 
consider as nucleated only those dEs having a 'significant' 
nucleus. Furthermore, and in order to minimise the impact 
of stellar mass effects, we only focus on galaxies in their 
'bright' dE subsample - roughly corresponding to systems 
having M t >2x 10 8 Mq . 

Table[T] summarises the main properties of metal-poor 
GCSs in these three dE subclasses, plus the field dlrrs of 
G10 with stellar masses larger than the already mentioned 
value. The non-diagonal elements in the upper part of the 
table contain the results of Kolmogorov-Smirnov tests for 
the corresponding T^^biue distributions, while the diagonal 
elements show the fraction of each population having ten 
or more blue GCs. The lower part of Table[T] provides a 
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Figure 2. Mass specific frequency of metal-poor GCs for cluster and field dlrrs, and Virgo dEs - including both nucleated and non- 
nucleated objects, as well as systems with disc features and central star formation. Galaxies having no blue GCs have been assigned a 
2jv blue = 0.12, while diagonal lines correspond to the loci of systems with 1, 10 and 100 blue GCs. The shaded region indicates, as in 
Fig.[TJ the potential location of harassed disc-like galaxies. Late-type field dwarfs preferentially have lower specific frequencies than both 
late- and early-type cluster dwarfs — a feature that is more prominent at the high-mass end and holds even for nucleated dlrrs. The 
arrow indicates the effect of M* doubling had the dE progenitors not had their star formation quenched (see text for details). 



more quantitative description of the distributions through 
their lower quartiles, their medians and their upper quar- 
tiles. The trend is unequivocal, and shows that the earlier 
the galaxy type, the richer its GCS. dE(N) contain, by far, 
the most populous systems and they appear to be inconsis- 
tent with those of field dlrrs at the 3 a level. The sole pres- 
ence of bright nuclei in dE(N) would already argue in favour 
of this distinction, but it is important to note that this differ- 
ence appears to hold for fainter nucleated dEs and dlrrs too 
(see Fig. [5| . dE(nN) and dE(di,bc) come next, and display 
rather similar frequencies. Table[T] shows that dE(nN), as 
also found by ML07, systematically have poorer GCSs than 
their nucleated counterparts - but their TN,blue distribution 
is still inconsistent with that of dlrrs at > 2 a level. The 
dE(di,bc) subclass is of particular relevance. They are char- 
acterised by the presence of spiral arms, bars and/or signs 
of recent star formation, and all these suggest they might 
be the most clear-cut examples of a recent environmental 
transformation. Indeed, Fig.[5]shows that a few of these sys- 



4 It is important to recall that the Kolmogorov-Smirnov test is 
specially sensitive to the median of the distributions, and hence 
the probability differences of these two subpopulations in Table[T] 



tems marginally occupy the region where we expect the ha- 
rassed remnants of more luminous disc-like galaxies to be 
found. Surprisingly, though, the rather elevated upper quar- 
tile of their frequency distribution, Qz(Tisi } blu&) = 18, indi- 
cates that a non-negligible fraction of these dEs have signif- 
icantly populous metal-poor GCSs. For instance, VCC 856 
-the paradigmatic example of a dE with spiral structure 
jjerien et al J l2000h - has Tn,mub = 16, a factor four larger 
frequency than the Sp with the highest Tjv,6i«e in Fig.[T] 

The paucity of metal-poor GCs in field late-type dwarfs 
represents a challenge for recent environmentally-driven evo- 
lutionary scenarios. Fading of infalling dlrrs after a recent 
(< 6 Gyr) gas stripping event does not change the number 
of GCs, and is thus incompatible with the rich systems we 
observe in M* >2x 10 8 M dEs - but appears to be fully 
consistent at lower masses. Could gravitational interactions 
then be responsible? We already pointed out that for harass- 
ment to increase the specific frequency of a given galaxy it 
is required that mass from the stellar body is preferentially 
lost with respect to the GC population. This can happen if, 
for example, GCSs have a similar or a more concentrated 
spatial distribution than stars do. This is a priori possible 
for low-mass late-type galaxies, as G10 show that the GCSs 
in dlrrs are essentially concentrated within the galaxy stel- 
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lar body (1 — 2r e ; see also lMora et al.ll2007T l. However, even 
if this leads to an increase of Tn,uub, harassment does not 
provide a satisfactory solution. It is clear from Fig. [5] and 
Table[T] that the number of metal-poor GCs in field dlrrs 
seems to be too low to reproduce the GCSs of most massive 
dEs even without invoking any kind of mass loss. The richer 
GCSs found in cluster dlrrs alleviate this tension, but imply 
that these two populations of star-forming dwarfs most likely 
had very different initial star formation histories - and that 
cluster dE p rogenitors are intrinsical ly different from field 
dwarfs (e.g., ISkillman fc Bendeilll995l ). 

4.2 GC Metallicity Distribution 

A second constraint can be derived from the GC metallicity 
distribution function of dEs and potential progenitors. It is 
well established that the GCSs of early-type galaxies present 
a bimodal metallicity (colour) distribution, with the frac- 
tion of metal-rich (red) clusters decreasing, and becoming al- 
most negligible, t owards fainter galaxy masses/luminosities 
|Peng et al.ll2006h . Disc-dominated L* galaxies also contain 
a relatively important metal-rich population that usually ac- 
counts for more than a third of their GCS. This is shown in 
Fig. [31 where we plot the fraction of red (metal-rich) clus- 
ters as a function of M* for the two samples (P08 and S08) 
having accurate estimates of subpopulations proportions 

The red GC fraction peaks at M + w 10 11 M , and de- 
creases towards both higher and lower masses. The majority 
(22/37) of dEs with M* < 2 x 10 9 M are consistent with 
having no red GCs at all. The wide range of metal-rich frac- 
tions displayed by dEs might suggest that these galaxies 
followed alternative evolutionary paths, or their progenitors 
were intrinsically different. 

4.2.1 LSBs vs dEs 

Therefore, while massive galaxies contain a significant popu- 
lation of metal-rich GCs, these red clusters appear to be less 
abundant in early-type dwarfs. In this context, for dEs to 
evolve from more massive systems it would require the latter 
to preferentially lose red GCs, so their relative fractions can 
be matched. This is in principle a difficult task as, at least in 
disc galaxies such as M31 and the Milky Way, the metal-poor 
popul ation is much m ore extended than the metal-rich one 
(e.g., iBica et al]|2006l ) - and LSBs appear to be consistent 
with this picture (V08). Given the strong galactocentric de- 
pendence of energy gain during an impulsive encounter, we 
would expect just the opposite trend -i.e., that the relative 
fraction of red clusters tends to increase. We again quan- 
tify this using the same MW scaling as in Sect. 14.170 at any 
given radius we estimate the fraction of enclosed disc mass, 
metal-rich ([Fe/H] > -1) and metal-poor ([Fe/H] ^ -1) 
clusters. 

The two blue curves in Fig. [3] show the evolution of 
the red GC fraction during this simple outside-in stripping 
model. The upper curve corresponds to an initial red GC 
proportion of ~ 30%, while for the lower one we randomly 
selected half of the original metal-rich subpopulation. It is 
clear that the more concentrated spatial distribution of red 
clusters leads to an enhanced relative fraction with respect 
to the reference value. If some dEs are the remnants of ha- 
rassed LSBs, they have to be found among the high red GC 
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Figure 3. Fraction of red (metal-rich) GCs as a function of host 
galaxy stellar mass. Galaxies consistent with having no red clus- 
ters have been slightly offset for displaying purposes. The gray 
curve is a running average, and shows that there is a trend for 
M* < 10 11 Mq galaxies to have lower metal-rich GC proportions. 
Blue curves indicate the expected evolution of a 10 10 Mq galaxy 
if mass stripping occurs from the outside-in (see text for details). 
In both cases the galaxy hosts a GCS with subpopulations having 
an identical spatial distribution to that of the MW, and initial 
red fractions of Ri30% and rj15%. Contrary to observations, the 
red GC fraction would tend to increase towards lower masses - 
their less extended spatial distribution makes them more resilient 
to stripping than the blue subpopulation. 



fraction specimens. The discrepancy might however be alle- 
viated if we consider that LSB progenitors can have a low 
bulge mass fraction - or be actually bulge-less. As metal-rich 
GCs are tradi tionally associated to the bulge component in 
disc galaxies (jForbes et al.|[200ll l. it would be possible that 
these systems possess a negligible fract ion of red GCs, as for 
example has been observed in M51 l)Chandar et al 1 12004J : 
but again see V08 for some evidence of red GCs in LSBs). 



4.2.2 dlrrs vs dEs 

A comparison with the GCS metallicity distribution of low- 
mass, late- type galaxies is perhaps more puzzling. G10 found 
a significant fraction of dlrrs with an abundant -and some- 
times even dominant- metal-rich subpopulation that could 
not be accounted for by contamination in their sample. 
Metal-rich GCs are not expected in low-mass systems from 
theoretical considerations on their formation process, and 
Fig. [3] indicates that they are actually not observed in a sig- 
nificant fraction of Virgo dEs. G10 suggest that very young 
ages or internal galactic extinction could be responsible for 
this effect. If metal-rich clusters turn out to be a significant 
subpopulation in dlrrs, it would be difficult to understand, in 
terms of environmental effects, how did these galaxies evolve 
into the metal-poor-dominated dEs we see in Virgo. Alterna- 
tively, the discrepancy between the metal-poor TN,biue val- 
ues of both galaxy types (Fig.0 might be somewhat allevi- 
ated if reddening is responsible for this effect. We however 
note that spectroscopic observations indicate that at least 
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several of these red candidates appear to be bona fide metal- 
rich GCs (I. Georgiev, private communication). 



4.3 Clustercentric Distribution 

We have so far focused our discussion on the GCS properties 
of Virgo dEs and their potential progenitors in the field. We 
now turn our attention to what can we learn from the char- 
acteristics of GCSs in early-type dwarfs within the cluster. 
Red circles in Fig.g] show T N>blU e for ACSVCS Virgo dEs 
as a function of their tridimensional distance to M87 (i.e., 
free from projection effects). As discussed by P08, there is a 
clear trend for specific frequency to decrease with radius out 
to ~ r2oo _ except for the two innermost, GC-poor dEs. One 
could argue that the trend is artificially created by contam- 
ination from the GC population of M87, but they show this 
is unlikely because i) the M87 contribution was accounted 
for when estimating the final GC specific frequencies, and 
ii) instead of being randomly distributed, GCs appear to be 
spatially concentrated around the galaxy body. 

In Sect . 13. f 1 we argued that harassment is an effective 
mechanism in a wide range of clustercentric regions, but 
that it is more efficient where the potential is deeper and 
in clusters with si gnificant substructure (e.g.. lGnedinll2003l ; 
ISmith et alJ l2010h . This is also indicated in Fig.0] where 
filled circles represent the fracti on of final stellar mass in 
the harassed dwarfs simulated bv lMastropietro et alJ (|2005T ) 
(their table f ) as a function of their time-averaged cluster ra- 
dius. As expected, the longer a galaxy resides deep within 
the cluster potential, the higher the mass loss is. In this con- 
text, the observed negative gradient in T^^iue towards the 
outskirts is puzzling and has in principle no trivial explana- 
tion within a framework of late accretion and recent evolu- 
tion through h arassment - as just the opposite trend is pre- 
dicted (see also lBekki et al.ll2003l ). If mass stripping is more 
efficient in the inner regions, how is it possible that precisely 
these dEs contain the richest GCSs? And what causes the 
GC deficiency towards the outskirts? One possibility would 
be that the outermost dEs have poorer GCSs because they 
have completed a larger number of orbits (were accreted 
earlier) within the cluster potential and are just caught at 
apoce ntre right now. This is however unlikely, as lGill et al.l 
(|2004 ) show that older subhaloes experience orbital circular- 
isation and preferentially end up closer to the cluster centre 
- and this appear s to be consistent with observations too 
jLisker et al J 12009 1. The two innermost dEs constitute the 
only exception to the radial trend, and, as discussed in P08, 
their GC paucity can be naturally explained within a grav- 
itational stripping scenario - be it through harassment or, 
most likely, due to the combined tidal fields of M87, dark 
matter and hot gas. 

A transformation via gas-stripping mechanisms seems 
even more problematic: why would progenitor dlrrs with 
more GCs per unit stellar mass have preferentially ended 
up close to the cluster centre? An intriguing possibility is 



5 (r) = a (1 — e 2 /2), where all orbits are assumed to be keplerian 
with semi-major axis a and eccentricity e. Note that the simulated 
dwarfs have been scaled so that the unstripped systems would 
appear to have Tjv blue. = 10, which roughly corresponds to the 
median value of actual Virgo dwarfs. 
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Figure 4. Empty circles show T N biue for ACSVCS Virgo dwarfs 
as a function of their tridimensional distance to M87, while the 
gray curve represents a smoothed version of the radial trend. 
Filled circles indicate the fin al fraction of stellar mass for the ha- 
rassed galaxies simulated by Mastropietr o et al . (2005) as a func- 
tion of their average clustercentric radius (see text for details). 
Galaxies that spend longer times deep within the cluster potential 
are affected the most. The trend of innermost dEs preferentially 
having richer GCSs than the outermost ones is the opposite than 
predicted by mass-stripping mechanisms. 



that a richer GCS is actually indicative of a higher DM-halo 
mass (Spi tler fc Forb es 2009), and we are just witnessing a 
total mass segregation effect. 



5 DISCUSSION AND FUTURE PROSPECTS 

In the previous Sections we have shown that GCSs can ef- 
ficiently constrain the role of the different physical mecha- 
nisms that operate in clusters and are supposed to shape 
early-type dwarfs in Virgo. The GCSs of low-mass (M* 
< 2 x 10 8 Mq) dEs indicate they are consistent with be- 
ing the remnants of gas-stripped late-type dwarfs, but leave 
very little room for mechanisms that involve any kind of 
stellar mass stripping. 

The origin of more massive dEs is however rather puz- 
zling. Their very high GC mass specific frequencies repre- 
sent a challenge for mechanisms that involve the removal of 
gas from infalling field dwarfs. Moreover, evolution through 
fast gravitational interactions from luminous progenitors ap- 
pears to be consistent with only a handful dEs exhibit- 
ing the lowest specific frequencies. Notably, one quarter of 
the dE(di,bc) subclass exhibit very high GC frequencies 
(Tjv,6iue > 18), a fact difficult to reconcile with the idea of 
a late transformation through harassment of accreted disc 
galaxies. This could be unex pected given prev ious studies 
suggesting the opposite (e.g.. iLisker et al1l2007f ). and these 
even include eviden ce from GCSs themselves. For instance, 
iBeaslev et al] (|2009h show that the GCSs of three Virgo dEs 
(VCC 1261, VCC 1087 and VCC 1528) exhibit significant ro- 
tational support, possibly indicating a discy origin for these 
systems. They propose that a combination of mass loss due 
to harassment and subsequent stellar population aging can 
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explain the high (v/cr) of their GCSs, and their compati- 
bility with the Tully-Fisher relation of late-type galaxies. 
Data from the ACSVCS however indicate that VCC 1261, 
VCC 1087 and VCC 1528 have T NMu e = 6, 19 and 23, re- 
spectively. These GC frequencies are 2-7 times larger than 
the typical values found in field disc galaxies (see Fig.[T|, 
so any mass loss scenario would require the progenitors to 
preferentially lose disc mass (with respect to GCs) by these 
factors in order to produce comparable remnants - an un- 
feasible scenario according to our analysis. 

The majority of these massive, GC-rich dEs are how- 
ever found within the nucleated subclass. These constitute 
an intriguing population with additional distinct properties, 
including the oldest stellar population ages, almost negli- 
gible rotational support, the roundest sha pes and a m ore 
concentrated clustercentric distribution (seejjLiskcr 200^ and 
references therein) . They comprise a significant fract ion of 
all dEs (~ 50-75%; ICote et alj|2005 iLisker et al.ll2007T ), and 
are traditionally thought to have formed along with the clus- 
ter. Having spent a considerable amount of time within the 
harsh cluster potential, it is remarkable that these systems 
have managed to retain such populous GCSs. This is most 
striking for the population residing in the inner cluster re- 
gion, where galaxies are supposed to move on ci rcular orbits 
suffer ing from substantial dynamical heating l|Lisker et al.l 
|2009f). The fact that dE(N) contain numerous and spatially 
extended GCs down to ~ 0.2 Mpc from M87 could be an in- 
dication that t hese systems are s hielded by a rather massive 
DM halo fe.g. JPennv et aDl2009l ) that has, so far, prevented 
the GCS disruption. In any case, their unique properties, in- 
cluding those of their GCSs, are totally inconsistent with an 
environmentally-driven late origin. 

A biased formation sc enario for Virgo dEs is however 
appealing in many aspects. iDiemand et al.l (|2005T ) show that 
the hierarchical assembly of DM haloes does not erase the 
memory of the initial conditions. As a result, the current spa- 
tial distribution and kinematics of surviving subhaloes es- 
sentially depend on the rarity of the primordial density fluc- 
tuations from which they formed. They find that rarer peaks 
-that depend on both the local and the underlying large- 
scale fluctuation amplitudes- collapse earlier and evolve to- 
wards a mo re strongly concentr ated configuration within the 
final halo. iMoore et al.l (|2006T l suggest that protogalaxies 
and metal-poor globular clusters form out of the first rarer 
peaks above a critical mass threshold for gas cooling. The 
resulting intense star formation activity led to reionisation, 
delaying the formation of subsequent baryonic structures - 
and suppressing it in the lowest mass haloes. Metal-poor 
GCs are believed to form in regions with high star formation 
densities and rates, and in direct prop ortion to the (ear ly) 
baryonic content of the (proto) galaxy (|McLaughlinlll 99*91 ) - 
wh ich is in turn relat ed to the (sub)halo mass. As discussed 
bv lPeng et alj (120081 ). in this scenario "GC formation is bi- 
ased toward the earliest collapsing halos that can create a 
large fraction of their stars at high efficiency before reioni- 
sation" . 

Remarkably, this picture provides a (qualitatively) ac- 
curate description of the bulk of observed properties of 
Virgo dEs. The existence of a m orphology-density rela- 
tion within dEs (|Lisker et alj I2007T ) is a natural outcome 
of this scenario. Indeed, the dE(di,bc)-dE(nN)-dE(N) se- 
quence can be thought of as originating from progressively 



rarer density peaks, thus explaining the more concentrated 
(extended) spatial distribution and prevalence of circular 
(radial) o rbits of nucleated (n o n-nucleated and d isc-like) 



dEs (e.g.. IConselice et alJliooH ; iLisker et all 120071 ). More- 



over, [Peni^^iL 1 20081 ) show that the negative radial gra- 



dient in Tjv for Virgo dEs can be qualitative ly explaine d 
in terms of a biased GC formation scenario l|Westfll993l ). 
Thus, the inner dEs would be systems that formed earlier, 
having higher star formation rates (SFR) and surface den- 
sities than dwarfs in the outskirts. This would in turn re- 
sult in a higher efficiency of GC formation, giving rise to 
the observed radial trend. Interestingly, high SFRs are also 
accompanied by energy release in terms of feedback from 
supernov ae and stellar winds, which c an create thick stellar 
systems llSanche z-Jansse n et al] 2010 |) and an extended GC 
distribution l|Mashchenko et all 120081 ), as observed in most 
Virgo dEs. Additionally, gas heating from these feedback 
mechanisms will make this component more susceptible to 
subsequent stripping even in the early stages of cluster for- 
matio n or during group pre-processing (e.g., ICortese et al.l 
2006). It is important to note that the majority of these 
dEs are too massive for feedback processes to completely 
expel their gas content, and environmental effects are there- 
fore needed to terminate the star formation activity and 
migrate the dwarf to the red sequence. In this sense, our 
analysis does not favour violent evolutionary mechanisms 
that result in substantial stellar mass losses, but more gentle 
processes involving gas remo val by a combinati on of inter- 
nal and external factors (e.g., Boselli et al. 2008) . The tida l 
stirring scenario l|Maver et al.ll200"ll ; lKazantzidis et al.ll201lh 
provides a plausible transformation mechanism, as long as 
GC and stellar mass losses are not prominent. 

The actual importance of stellar mass loss in cluster 
galaxy evolution can be quantified through the study of 
intra-cluster globular clusters (IGCs). It is now well estab- 
lished that a significant population of such unbound systems 
inhab it the intergalactic regions of re l atively massive c lus- 
ters dLee et al.ll201(il ; iPeng et al.ll201ll ; IWest et al.ll201ll ). It 
is however not yet clear what fraction of these IGCs originate 
from tidal interactions and what fraction co nstitute a pri- 
mordial cluster population l|West e t al. 1995). In particular, 
spectroscopic determinations of velocities and metallicities 
for these systems will prove critical to address their dynam- 
ical s tatus and their connectio n with cluster galaxies GCSs 
(e.g.. [Romanowskv et al.ll2012l ). 

We finally stress that the analysis of GCSs suggests 
that the progenitors of Virgo dEs are intrinsically differ- 
ent from today's field dwarf galaxies - but cluster dlrrs ap- 
pear to be fully consistent. Interestingly, the significantly 
different metal-poor GC frequencies of these two late-type 
dwarfs can also be naturally explained within the biased for- 
mation scenario. Virgo dlrr galaxies collapsed earlier from 
higher-CT peaks than their counterparts in low-density en- 
vironments, thus forming a larger fraction of their current 
stellar mass in the form of GCs before reionisation. Star for- 
mation was resumed shortly after this event, but under very 



6 The sole existence of the star-forming dE(bc) subclass rules 
out the former possibility, while the relevant role of environment 
is reinforce d by the existence of a morphology-density relation 
within dEs iLisker et al.ll2007ll . 
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different conditions for both systems: while cluster dlrrs had 
to struggle within the developing harsh potential well, field 
dwarfs have been able to form stars without significant dis- 
tress until today. The difference in Tjv.Wtte further suggests 
that, modulo internal evolutionary effects in the GCS, the 
initial star formation histories of cluster and field dwarfs 
m ust also differ. This nothin g but reinforces the conclusion 
of lSkillman fe Benderl | | 1995ft that "present day dEs and dls 
may share a common ancestor, but dEs do not evolve from 
dls". 

In this study we have shown that GCSs are very sen- 
sitive tracers of evolutionary mechanisms in high-density 
environments, and thus put strong constraints on the rel- 
evance of the various processes. In forthcoming papers we 
will present a more elaborated treatment of the effects of 
fast interactions on the GCSs of early- and late-type galax- 
ies using numerical simulations, which will further allow us 
to compare their kinematics with available observations. De- 
spite of the growing amount of data, our understanding of 
GCSs in Virgo dEs is still rather limited. The situation will 
dramatically change once the Next Generation Virgo Survey 
(Ferrarese et al. 2012) is completed. The analysis of GCSs 
in a complete sample of Virgo dwarfs as a function of their 
stellar mass, morphological type and clustercentric position 
should provide us with a clear picture of their origin and the 
role of environmental effects in shaping their properties. 
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APPENDIX A: A NOTE ON THE ORIGIN OF 
CLUSTER SO GALAXIES 

Figures[T] and [3] show that SO and Sp galaxies share similar 
locations in both the M + -Tjv,biue and M*-red GC fraction 
planes. Taken at face value, this would argue in favour of 
an evolutionary link between both galaxy types, as envi- 
ronmental effects probably do not change significantly the 
quantities involved: ra m-pressure strippin g does not affect 
the stellar components jQuilis et alj 1200(1 ) and harassment 
leads to very little disc mass loss for these massive systems 

jMoore et alJll999h. 

lAragon-Salamanca et all (|2006l ) were the first to pro- 
vide a quantitative comparison between SOs and Sps based 
on the properties of their GCSs. They find that Sn is on 
average a factor ~ 3 larger for the former population than 
for the latter, and suggest that, under the premise that the 
number of GCs does not vary, it can be interpreted as ev- 
idence for luminosity fading due to star formation shutoff 
in SOs. The constant N gc hypothesis between the two pop- 
ulations is however questionable, as at least 2/8 of the SOs 
in their table 1 have more GCs than any of the Sps - and 
the discrepancy is even larger if Sombrero (NGC 4594) is not 
considered to be a spira l (e.g., Gadotti fc Sanchez- Jansse"r] 
120111 ; ISpitler et all l2008i : iRhode fc Zepj|2004l ). The fading 
interpretation might thus not be that straightforward for 
the most massive SOs harbouring populous GCSs - but is in 
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principle consistent with the origin of less luminous lentic- 
ulars. This further supports the existence of two different 
formation mechani sms for central and sate llite SOs, as re- 
cently proposed bv lWilman fc Erwinl (|2012h . 

While a detailed comparison between the GCS proper- 
ties of SOs and Sps is beyond the scope of this paper, we note 
that Fig.[T] also shows a slight hint of SOs having a higher 
mass specific frequency than late-type disc galaxies. Uncer- 
tainties in both stellar masses and the number of GCs are 
however significant compared to the differences, and a more 
detailed study complemented with independent observables 
(e.g., bulge-to-disc mass fractions, stellar population ages 
and metallicities) is highly desirable in order to clarify the 
origin of SOs. 
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